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Abstract 

Electrodynamic spherical harmonic is a second rank tensor in three- 
dimensional space. It allows to separate the radial and angle variables 
in vector solutions of Maxwell's equations. Using the orthonormaliza- 
tion for electrodynamic spherical harmonic, a boundary problem on a 
sphere can be easily solved. 

1 Introduction 

In this paper we introduce new function — electrodynamic spherical har- 
monic. It is represented as a second rank tensor in three-dimensional space. 
But the function differs from tensor spherical harmonic [1—5]. By its defini- 
tion the electrodynamic spherical function possesses a number of properties 
of usual scalar and vector harmonics and includes them as component parts. 
Why electrodynamic spherical harmonic? The word "electrodynamic" im- 
plies that the function is applied for solution of vector field problems. We 
use it in electrodynamics, however one can apply the function for description 
of spin fields in quantum field theory. 

Electrodynamic spherical harmonic is not a simple designation of the well- 
known functions. It satisfies the Maxwell equations and describes the angular 
dependence of vector fields. The introduced function separates the variables 
(radial coordinate and angles) in the fields. Moreover, the notation of the 
fields in terms of electrodynamic spherical harmonic noticeably simplifies the 
solution of a boundary problem on spherical interface. Just application of 
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the ortonormalization condition allows to find the coefficients of spherical 
function expansion (e.g. scattering field amplitudes). 



2 Scalar and vector spherical harmonics 

Hamilton's operator V in three-dimensional space contains the derivatives on 
three coordinates. In spherical coordinates (r, 9, tp) the derivatives on radial 
coordinate and angles can be separated. This is achieved by representing the 
unit tensor 1 in 3D space as the superposition of two projection operators: 

(r 6?) r r x 2 \ 1 1 

^-7rJV = ^r(rV)--rX(rXV), 

where r ® r/r 2 is the projector onto the direction e r = r/r, — r x2 /r 2 is the 
projector onto the plane orthogonal to the unit vector e r . Tensor r x dual to 
the vector r gives the well-known vector product when acting on a vector a: 
r x a = r x a and ar x = a x r [6]. Introducing the vector differential operator 

L = yr x V (1) 



the equation (??) is rewritten as follows 



V = -#- - \r x L. (2) 

r or r z 

Vector L is called orbital angular momentum operator in quantum me- 
chanics, because it is presented as vector product of radius vector r and 
momentum p = — iV operators. L includes only derivatives on the angles 9 
and ip. Using the definition of L one can derive the Laplace operator 

(3) 

/ Ul \ Ul I I 

and the following properties: 

rL = 0, (Lr) = 0, L 2 L = LL 2 , L x L = iL. (4) 

Scalar spherical harmonic Yi m (9, <p) is defined as the eigenfunction of the 
operator L 2 : 

L 2 Y lm = 1(1 + l)Y lm , (5) 
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where I and m are integer numbers. Number m is the eigenvalue of the 
operator of projection of angular momentum onto the axis z, L z : 

L z Yi m = mY lm . (6) 

Spherical harmonics Y\ m are orthogonal and normalized by the unit: 

r r Y Vm>& <p)Yi m (0, <f) sin 6d6dip = 8 vl 5 m , m . (7) 
Jo Jo 



where the sign * denotes the complex conjugate. 

If we multiply equation (j3J) by the vector operator L and take into account 
the commutation of L and L 2 , we will obtain that vector ~LYi m satisfies the 
same equation (j5J), too. The quantity defined as 



1(1+1 



=LYz m (8) 



is called vector spherical harmonic. The coefficient before LY/ m is chosen so 
that the orthonormalization condition is of the form 

/■7T /*2tT 

/ / X*, m ,(0, tp)Xi m (8,<p) sm6d6d<p = 5i>i8 m >m- (9) 
Jo Jo 

From the self-conjugacy of the angular momentum operator L and prop- 
erties (0J the orthogonality 

f I** e r {lC Vm , x X, m )sin0d0d^ = f f^^" 1 ^ sinfldfldy 
Jo ' Jo Jo ^(z + !)//(// + i) 

7 ° 7 ° v /^ + l)z / (£' + 1) 

follows. Below we give some properties of vector spherical harmonics: 



LX im = + l)Y lm , L(e r xX /m ) = 0. (11) 

Scalar (vector) spherical harmonics satisfy the scalar (vector) equation for 
eigenfunctions of the squared orbital angular momentum operator L 2 . The 
orthogonality condition for the scalar (JTj) and vector (jUJ) spherical harmonics 
have the same form. Therefore, one can hope to combine them into one 
mathematical object. 
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3 Electrodynamic spherical harmonic: defi- 
nition and properties 

We define an electrodynamic spherical harmonic as a second rank tensor in 
three-dimensional space 



The first term of (IT2|) determines the longitudinal part of the tensor. It is 
calculated by means of the scalar spherical function. The last two summands 
of (fT2|) fix the transverse solution, in the plane (8, if) perpendicular to the 
direction e r . It includes vector spherical harmonics. The left and right 
vectors in dyads form two sets of orthogonal vectors: (F; m e r , X im , e r x X im ) 
and (e r ,e ,e^). 

3.1 Orthonormalization 

Multiplying the electrodynamic spherical harmonic by the Hermitian conju- 
gate we get to 

Fih'Ftm = V m /Vj m e r (S)e r + X*, m ,Xj m (e e (8)e 6 ) + e ¥ ,(8)e ¥ ,) + (e r (X i * m , xX lra ))e r x . 



The quantity before each dyad is orthogonal or normalized as ©, (J9]), or 
({TO]) . Therefore, the orthonormalization condition for electrodynamic spher- 
ical harmonics is 



Here we consider the vectors e r , e e , and e 9 to be constant in dyads and dual 
tensor e* . If tensors F[ m are the solutions of equations, we can always write 
these solutions in components. For each component of the tensor F^F^ the 
orthonormalization is carried out. We will obtain the same, if the dependence 
of the orts e r , eg, e v on angles in (fT3|) is omitted, i.e. they are regarded as 
constants. However, when substituting F [m in equations, we should take into 
account the angular dependence of the orts. 



Fim = Y lm e r ®e T + X lm <g> e e + (e r x X lm ) <g> e, 



(12) 



(13) 




(14) 
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3.2 Explicit form 

Let us substitute the explicit expression for the vector operator L 



d 



-ie„ 



e e d 



(15) 



""89 ' 'sinddtp 

into equation (IT2"|) . Then the electrodynamic spherical harmonic is equal to 



lm 



e r ® e r + 



sm9Jl(l + l) V d( P de , 



Y 



Imi 



(16) 



where / = — e*e f x = 1 — e r <g> e r is the projection operator onto the plane (9, 
ip). To calculate the derivatives one can replace them by means of operators 
L z and L± = L x ± L y as 



e iv L_ 



because their action on the scalar spherical harmonic is well-known: 

L + Y im = yj(l-m)(l + m + V)Y hm+x , L_Y lm = y/(l + m)(l-m + l)lj, m _i. 

(17) 

Hence, the electrodynamic spherical harmonic can be presented as follows 



lm 
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e r (g> e r 



sin ej/(z + 1; 



iX* + e: 



i sin 6 



-up ; 



e iv L_ 



(18) 

It is easy to exclude the operators from (jTSJ) . The final formula for tensor 
F; m contains scalar spherical harmonics as angle dependence. Unit vectors 
e r , eg, determine the structure of the tensor in three-dimensional space. 
Fi m is formed by three basic tensors: e r <g> e r , e*, and /. Hence, it commutes 
with each of these tensors. 



3.3 Invariants 

The first invariant of the electrodynamic spherical harmonic as tensor quan- 
tity is its trace. The trace of the tensor (1T21) equals 

tr(F lm ) = Y lm + 2(e e X lm ). (19) 
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The second invariant is determinant 

det(F lm ) = Y lm X.* m . (20) 

The third invariant of three-dimensional tensor Fi m is the trace of the ad- 
joint tensor F im . Adjoint tensor is defined by F im Fi m = Fi m Fi m = det(F im )l 
and equals 

Fi m = Xf m e r (8 e r + Yi m e e <g> Xi m + Yi m e v <8> e r x X !m . (21) 

Further the trace is easily calculated: 

tr(F, m ) = XL + 2F /m (e e X lm ). (22) 

Using these three invariants one can find other ones. For example, the trace of 
squared tensor is determined from equation tr(F^) = (tr(F; m )) 2 — 2tr(F/ m ). 

3.4 Generalization of electro dynamic spherical harmonic 

The main condition on electrodynamic spherical harmonic is the orthonor- 
malization (JHj). There is more general form of the second rank tensor spher- 
ical harmonic satisfying this equation. It is 

Gim = Y lm e r ® a + X im <g> b + (e r x X im ) ® c. (23) 

Unit vectors a, b, and c form the orthogonal basis in three-dimensional space. 
In equation ( 1121) we have assumed a = e r , b = e$, c = e v because of the 
spherical symmetry of the function. If we will take a = e r , b = e#, c = — e v , 
then the electrodynamic spherical harmonic becomes more complex function 
in explicit form, however its invariants are simplified (e.g., the trace equals 
tr(F lm ) = Y lm ). 

4 Solution of Maxwell's equations 

In this section we will find the spherically symmetric solutions (i.e. electric 
E and magnetic H fields) of the Maxwell equations 

V x E(r) = ifc/iH(r), V x H(r) = -ifeE(r) (24) 
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for the monochromatic electromagnetic waves in isotropic medium with di- 
electric permittivity e and magnetic permeability \i. The quantity k = u/c is 
called wavenumber in vacuum, and u is the wave frequency. Arbitrary time 
dependence can be obtained by using the linear superposition of monochro- 
matic waves E(r) exp(— iuot). 

We will search the solution in the form 

E(r) = F lm (0, V )E l (r). (25) 

The components of the vector E ; in spherical coordinates depend only on the 
radial coordinate r. The dependence on the angle coordinates presents only 
in the basis vectors. So, the vector E' looks like 

E'(r) = E l r (r)e r + E l e (r)e e + £j,(r) V (26) 

Further we should calculate rotE. By substituting Hamilton's operator 
(EJ) one obtains 

<9E i i i 
V X E = °*dF " r L(GrE) + r Gr(LE) ' (27) 

where the arrow j implies that operator L acts only on vector E, but not e r . 
Let us calculate each summand of equation (1271) applying the solution (J2"5]l . 
The first term is of the form 

«ff = F <~<TF' (28) 

where the commutation relation [Fi m ,e*] = is taken into account. The 
second summand yields 

L(e r E) = L(e r E) - L(e r E) = L(Y lm E l r ) - L(e r E). (29) 

The quantity L(e r E) is easily calculated using the explicit expression (1T51) of 
the operator L and the relationships de r /d8 = eg and de r /dip = e^sin^: 



L(e r E) = -ie T x E. (30) 



So, we get the formula 
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L(e r E) = F lm J 1(1 + l)e fl ® e r + ie r x E'. (31) 



The third term in (j27p can be rewritten using the equation (1111) : 



e r (LE) = e r (E l r L(Y lm e r ) + E l g LX lm + ^L(e r x X lm )) = y/l(l + l)F lm (e r ®e e 

(32) 

Finally, the curl of the electric field vector E equals 



d_ , 1 y W + l) 
dr 



V x E = F lm (9,<p) I e?^ + -e r >< - v \ e* ] E'(r). (33) 



In expression (1331) we have took into account the derivatives on the angle 
variables. Therefore, further the orts of spherical coordinates e r , eg, and 
can be considered as constants. The Maxwell equations are reduced to the 
set of ordinary differential equations 



dE l 1 , + 1) 



e 



? jg + l e >H'- i V' ( ' + 1) e-H- = - ite E'. (34) 
dr r r r 



Equations ( jS|) allow to determine the radial dependence of the fields. 
Multiplying the set of equations AM]) by the unit vector e r we can express 
the longitudinal components of the fields as follows 



^ " — ^' ^ " ~Hh— Ee - (35) 

The tangential field components E{. = IE' = E l e e e + E l and H' = IH l are 
determined from the equations which can be written in matrix form: 

d(rW) 

— - = ikMirW), (36) 
dr 

where 

W-(3). -=(-l A-^-g^.^ 

Equation fl36|) is satisfied for inhomogeneous media e(r), /■*(?")> too. Such 
matrix equation can be solved numerically for arbitrary medium, or analyti- 
cally for homogeneous one. Let us find tangential field components W when 
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e = const and \i = const. The simplest way is to write the equation for 
projection W e = (H l 9 ,E l 9 ): 



dr 2 \ r 



The solutions of the equation (1381) are well-known and can be presented as 

W e = /« ( * ) + ^ = (/ (1) Cl + f^c 2 ) ( * ) , (39) 

where Ci and c 2 are constant vectors. The couples of independent solu- 
tions are spherical Bessel functions = ji(nkr), = j_i_i(nkr) or 
spherical Hankel functions of the first and second kind /W = h^\nkr), 
/(2) = hf\nkr). n = y/ejl is the refractive index. After determining the 
99-projections of the fields as 

Wip kre/j ( li ) dr 
one can write the transverse vector field 



where 



/ifcr dr 



771,2 = / (1,2) eg <8> e e - - 7 - — <g> e^, 



0,2 = / (1 ' 2) e e ® + —^4 ^ ® e». (42) 

ekr dr 

Tangential field vector W plays an important part, because it is contin- 
uous on the spherical surface. That is why it can be applied for the study of 
electromagnetic wave diffraction by a sphere. 



5 Conclusion 

Thus, the general solution of the Maxwell equations is of the form 

( H(r) U « £ Flm{e ^v\r) ( #| )( * V (43) 
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where V is the matrix that takes into account the longitudinal components 
of electric and magnetic fields. This matrix can be easily calculated from 
the equation fl35l) . In each partial solution included into the general one f|43l) 
the radial and angle variables are separated. Using the orthonormalization 
f fT4|) for the electrodynamic spherical harmonic Fi m , each partial wave can 
be easily singled out: 

v i M ( V[(r) vi(r) )( c\\ _ r + ( H(r,0'V) \ q'aq'a.J 

V (0 { C[(r) (l(r) )[ci)-J Jo m ' * } [ E(r, 9', ) ^ 68 " 

(44) 

This property of the electrodynamic spherical harmonic is very useful for the 
investigation of electromagnetic wave scattering. In scattering, the boundary 
condition is the single equation for tangential fields W. It is easily solved, 
if the orthonormalization is applied. Some attempts of investigation of scat- 
tering in the similar manner as described above have been made in [7]. 

In the general solution (|4~3I) the constants c 12 determined by initial condi- 
tions are explicitly shown. Vectors Ci and c 2 set independent solutions. For 
instance, if Ci = 0, then the radial dependence is determined by the function 
f( 2 \r), and vice versa. If c 2 = and f^(r) = h^(nkr), then the field f|43l) 
is the multipole expansion [8]. The amplitude of electric multipole field is 
equal to a£;(/,m) = c^, the amplitude of magnetic multipole field is equal 
to a^j(l, m) = Cie v . So, vector ci can be called vector amplitude of multipole 
fields. 

In further investigations we will study the scattering and multipole ex- 
pansion of electromagnetic fields in details. 
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